rm(list = ls())

load("st_mn0_100.RData")
load("gt_mn0_100.RData")

c1 <- apply(res_gt_mn0_100, 1, function(i) sum(i=="reject"))/1000
c2 <- apply(res_st_mn0_100, 1, function(i) sum(i=="reject"))/1000

load("st_mn9_100.RData")
load("gt_mn9_100.RData")

c3 <- apply(res_gt_mn9_100, 1, function(i) sum(i=="reject"))/1000
c4 <- apply(res_st_mn9_100, 1, function(i) sum(i=="reject"))/1000


load("st_mn0_150.RData")
load("gt_mn0_150.RData")

c5 <- apply(res_gt_mn0_150, 1, function(i) sum(i=="reject"))/1000
c6 <- apply(res_st_mn0_150, 1, function(i) sum(i=="reject"))/1000

load("st_mn9_150.RData")
load("gt_mn9_150.RData")

c7 <- apply(res_gt_mn9_150, 1, function(i) sum(i=="reject"))/1000
c8 <- apply(res_st_mn9_150, 1, function(i) sum(i=="reject"))/1000

cbind(c1,c2,c3,c4,c5,c6,c7,c8)

